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Abstract 

In a dedicated test setup at the Kamioka Observatory we studied pulse shape dis- 
crimination (PSD) in liquid xenon (LXe) for dark matter searches. PSD in LXe was 
based on the observation that scintillation light from electron events was emitted 
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over a longer period of time than that of nuclear recoil events, and our method used 
a simple ratio of early to total scintillation light emission in a single scintillation 
event. Requiring an efficiency of 50% for nuclear recoil retention we reduced the 
electron background to 7.7±l.l(stat)±Q;|(sys) x 10 -2 at energies between 4.8 and 
7.2keV ee and to 7.7±2.8(stat)±| j(sys) x 10~ 3 at energies between 9.6 and 12keV ee 
for a scintillation light yield of 20.9 p.e./keV. Further study was done by masking 
some of that light to reduce this yield to 4.6 p.e./keV, the same method results in an 
electron event reduction of 2.4±0.2(stat)±Q2(sys)xlO _1 for the lower of the energy 
regions above. We also observe that in contrast to nuclear recoils the fluctuations in 
our early to total ratio for electron events are larger than expected from statistical 
fluctuations. 

Key words: Scintillation, Liquid xenon, Pulse shape discrimination 
PACS: 29.40.Mc, 23.40.-s 



1 Introduction 



The results of various astronomical observations [T] - [5] show strong evidence 
for a large amount of dark matter in the universe. Weakly Interactive Massive 
Particles (WIMPs) are a dark matter candidate motivated in extensions of the 
theory of high energy particle physics [6]. Various dedicated WIMP dark mat- 
ter search experiments are underway around the world [7] - [9]. The XMASS 
experiment, using liquid xenon (LXe) as a target for WIMP dark matter, was 
proposed in 2000 [10] . The construction of the 800 kg detector was finished in 
2010. 

The interaction of WIMP dark matter is observed as a nuclear recoil in a 
specific detector's target material, which in our case is LXe. The main back- 
grounds (BG) for such nuclear recoil events are electron events (from pho- 
toabsorption or Compton scattering of environmental gamma rays), nuclear 
recoils from the scattering of fast neutrons, and possibly alpha and beta de- 
cays in the detector medium itself. The aim of this study is to use the shape 
of scintillation light pulse in LXe to discriminate against electron events. 

The scintillation mechanism is classified into two processes determined by 
whether or not the process includes recombination [II] . 

the process without recombination 
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Xe* + Xe^Xe*. 

Xe*. -> 2 Xe + hu 



(1) 



the process with recombination 



Xe+ + Xe^Xe^ 
Xe+ + e"^Xe** + Xe 
Xe**^Xe* + heat 
Xe* + Xe^Xe* ; 

Xe*. -> 2 Xe + /«/ 



(2) 



Both processes lead to the formation of exciton ( Xe*. ). The de-excitation has 
two components called singlet and triplet component. The singlet component 
is caused by a spin singlet state( 1 £+ ), and the triplet component is caused 
by a spin triplet state( 3 S„ ) [T2"] . 

The shape of the scintillation light pulse in LXe is determined by the lifetimes 
of the excited states, the time scale of electron-ion recombination, the timing 
resolution of the photo detector, the time of flight in the detector, and elec- 
tronics employed to record the scintillation light. The convolution of all these 
components shapes the recorded scintillation light pulses. 

Pulse shape measurements for 1 MeV electrons, a particles and relativistic 
ions in LXe were reported [13]. In the case of a particles, two distinct lifetime 
components were observed. The lifetimes for singlet and triplet states were 
found to be approximately 4 ns and 22 ns respectively. In the case of electron 
events on the other hand, only one component with 45 ns lifetime was observed. 

This pulse shape difference is attributed to the influence of electron-ion re- 
combination [TJ] [15] . In the case of electron event, dE/dx is one order of 
magnitude smaller than that in the case of a nuclear recoil [H] [16], and the 
ensuing ionization is thus spread out over a larger volume. Therefore the re- 
combination process of electron events takes longer than that of nuclear recoil 
events and dominates the pulse evolution. This hypothesis was confirmed by 
measurements in electric fields, as the slow component was not observed when 
an external electric field is applied. Under the influence of an external electric 
field recombination is suppressed by the drifting apart of the opposing charges, 
and the two components characteristic of high dE/dx events re-emerge, so that 
the pulse shape of electron events in an electric field of 4kV/cm is very close 
to that of nuclear recoils in zero electric field [12]. In this paper we aim to 
exploit the characteristically long time constant observed for electron events 
in LXe at zero electric field. 
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To this end we examine pulse shape discrimination (PSD) in LXe at energies 
of less than 20keV ee ( electron equivalent keV). This energy range is most 
relevant for dark matter searches. The average waveform of nuclear recoil and 
electron events was previously compared above 10keV ee [T7]. But no event- 
by-event analysis was done. An event-by-event analysis was reported [18J, but 
the light yield was one order of magnitude lower than our setup. Using dual 
phase detectors the PSD study was also performed [19] . A reduction factor 
for electron events of 0.2 by PSD alone was demonstrated at 5keV ee in a dual 
phase detector by limiting the electric field to 0.06kV/cm. 

Previous investigations ([T7] - [TO]) were done at scintillation light yields be- 
low 5 p.e./keV (p.e. photoelectrons) . Using two closely spaced photomultipliers 
(PMTs) of the kind also used in the XMASS experiment we measure a light 
yield of 20.9 p.e. /keV, which is expected to us to clearly observe a difference 
in the pulse shapes of nuclear and electron events even at energies as low as 
5keV ee . Anticipating more limited light collection in real detectors we artifi- 
cially reduced the photosensitive area in our experimental setup and repeated 
the measurements at an effective light yield of 4.6 p.e./keV. 



2 Detector setup 

Our measurements were made in a dedicated setup shown in FigJT] Two 2 inch 
hexagonal PMTs (Hamamatsu R10789) facing each other from a distance 
of 6 cm are viewing the intervening LXe volume from the top and bottom 
respectively. To the sides the volume is limited by a highly reflective PTFE 
surface enclosing 0.58 kg of LXe in total. Embedded in the PTFE reflector 
were a LED and a 300 Bq 57 Co source. The LED was used to obtain single 
p.e. spectra and the source to monitor the PMT gain. The light yield from this 
source was found to be stable within 2% and is incorporated in our systematic 
uncertainty. To protect the LXe from radioactive contamination this 57 Co 
source was enclosed in a thin stainless steel container. 

PMTs, LXe, and the calibration light sources described above were all kept at 
-100°C inside a stainless steel vessel that itself is suspended inside a vacuum 
chamber for thermal insulation. Radioactive sources placed outside of the LXe 
volume and its vacuum enclosure produce the recoil events used in this study: 
Either a 137 Cs source provides gamma rays that produce electron events or a 
gamma tagged 252 Cf fission source provides neutrons to study nuclear recoils. 
To tag the fission events a plastic scintillator and a PMT (PMT3) were set 
up next to the Cf source. In the direction of the LXe on the other hand the 
gamma rays emitted in the 252 Cf fission events were shielded by 5 cm of lead. 

When filling the setup with LXe the Xe gas was passed through a SAES 



3 



getter (Model PS4-MT3-R-1) to remove impurities in the Xe gas. Beyond 
this no further efforts were made to clean up the high purity Xe used in this 
experiment. 

For data acquisition we used NIM logic to trigger recording of the waveform 
in a LeCroy WavePro 900 digital oscilloscope as shown in FigfJJ Events in the 
LXe volume are identified by a coincidence of signals in the two PMTs that 
view the LXe volume. The width of the coincidence timing window is 100 ns, 
and the discriminator thresholds were set to 2p.e. equivalent. As mentioned 
above a plastic scintillator is used to tag 252 Cf fission events every time the Cf 
source is employed to produce neutron events in the LXe. The corresponding 
signal from the PMT that reads out the plastic scintillator starts a 140 ns gate 
that is only for neutron recoil data is also entered into the coincidence unit. 

The digital oscilloscope provides 8 bit resolution and a 1 GHz sampling rate. 
We used it to record 10 /j,s traces of all PMT signals involved in the respective 
measurement with 1 ns timing resolution. 

The gain of the two PMTs reading out the LXe volume were both set to 
3.8xl0 6 at the operating temperature of -100 °C, with the HV for PMT1 and 
PMT2 fixed at 1.26 kV and 1.21 kV, respectively. Using the lp.e. calibration 
obtained from the LED data, the fit shown in Figj3]to the 57 Co spectrum deter- 
mines the energy resolution to be 5.4% (RMS) at 122 keV and the scintillation 
light yield to be 20.9 p.e./keV ee . 

Light collection in real dark matter detectors will be more limited than in 
this small scale dedicated setup. To quantify the impact of the ensuing loss 
in statistical power on our discrimination method we artificially reduced the 
light collection in our setup by covering part of the photocathode area of the 
bottom PMT (PMT2) with a copper mask. While the PMT1 signal was no 
longer used to estimate the deposited energy in our study of reduced effective 
light yield, its signal was used in the trigger during that study in just same 
way as usual. 

The resulting 57 Co spectrum is shown in FigJU using the p.e. count of only 
this bottom PMT. The effective light yield was thus reduced to 4.6p.e./keV ee 
and the energy resolution to 11.7% (RMS) at 122 keV. 



3 Data reduction 

A trigger offset of 260 ns allowed us to monitor the baselines 260 ns prior to 
the recording of our physics events. To clean up both our electron and nuclear 
recoil event samples we require that none of the events show any "pre-activity" 
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in the first 150 ns from the beginning of its PMT traces, i.e. from 260 ns to 
110 ns prior to the event trigger time. 

The dynamic range of oscilloscope was chosen to saturate at a signal height 
of -400 mV for the PMTs that read out our LXe volume. To avoid problems 
with saturation, we disregard events where PMT traces from either PMT1 or 
PMT2 surpass 350 mV, which is roughly equivalent to 1000 p. e. for electron 
events. 

Neutron events are identified by the neutron's time of flight (TOF) in the 
offline analysis. Fig[5] shows the TOF distribution extracted from software 
threshold crossings between the leading edges of PMT3 and PMT2. The soft- 
ware thresholds for all timing determinations corresponds to the typical height 
of a 3p.e. signal from the PMTs. The peak at zero timing is due to 252 Cf fis- 
sion gammas passing through the lead shielding set up to suppress such direct 
gamma interference. To limit the gamma ray background in the sample and 
specifically select fast neutrons a narrow TOF range from 15 ns to 30 ns was 
chosen for further analysis. 

Table [TJ summarizes our data reduction: 





252 Cf (neutron) Run 


Cs (gamma ray) Run 


before cuts 


8.0xl0 5 


4.0xl0 5 


after pre-activity cut 


5.6xl0 5 (69%) 


3.5 xlO 5 (88%) 


after saturation cut 


2.2 xlO 5 (27%) 


7.4xl0 4 (18%) 


after TOF selection 


2.2 xlO 4 (2.7%) 





Table 1 

Data reduction summary of PSD measurement. 



4 Pulse shape discrimination 

The parameter we chose to discriminate between electron and nuclear re- 
coil events is the ratio Rpsd of scintillation light detected in just the first 
At t i =20 ns of the scintillation light pulse, which we refer to as the prompt 
light, to the total amount of scintillation light detected in that same pulse. 
As we measure light with PMTs, the amounts of prompt and total light are 
recorded in units of p.e.: p.e. prmt and p.e. tat . The window for evaluating the 
amount of total light starts at the same time to as that for the prompt light 
and is At to t =200 ns long. t is the same software threshold crossing time as 
used in the TOF distribution above. Therefore each PMT has its own soft- 
ware threshold crossing resulting in potentially slightly different start times 
to,i and to,2 for the integration of a PMT's respective scintillation signal. The 
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respective integrals for prompt and total charges measured from the base line 
subtracted oscilloscope traces VpMTiit) and VpMT2{t) of the scintillation signal 
recorded by PMT1 and PMT2 respectively are added as the ratio is evaluated 
for each event: 



R Jtp.i VPMT1 dt + J to ° 2 2 " VpMT2 dt _ p.e.prmt , . 

Jto ! VpMTl Ot + J t 2 t / PMT2 «C 



The 20 ns width for the prompt timing window was optimized for best dis- 
crimination against gamma recoils. In the case of the reduced effective light 
yield measurement all terms involving Vpmt2 are simply dropped. A possible 
systematic effect of the recoil events position in the LXe volume was studied 
with the help of the light balance in the two PMTs. 



4-1 Neutron data 



The recoil energy E r of a xenon nucleus elastically scattered on by a neutron 
with energy E n is expressed by the following equation: 



2(A+1- cos 2 9 - cos 9VA 2 - 1 + cos 2 9) 

Er = En {iTat (4) 

where A is the mass number of the recoil nucleus, and 9 is the scattering angle 
of the neutron. The maximum recoil energy at 9 ~ 180° is 220keV r for an 
8 MeV neutron. To estimate the light yield we need to know the relative scin- 
tillation efficiency L e ff that describes the scintillation light yield of a nuclear 
recoil event as compared to the yield of an electron event at the same energy; 
we assume L e ff = 0.2 [9]. With that the maximum visible energy for a Xe 
nucleus recoiling in LXe from the elastic scattering of an 8 MeV neutron is 
44keV ee , corresponding to about 920 p. e. in our setup. 

The trigger rate for the neutron data run with the 252 Cf source was 13.2 Hz (co- 
incidence rate of PMT1&PMT2&PMT3). The trigger rate of PMT3 alone was 
16.2 kHz and the coincidence of just PMT1 and PMT2 occurred at 1.12 kHz. 

Accidental coincidences in the sample can be estimated from the background 
in the TOF distributions for each energy range individually, and are found to 
be less than 5% after our event selection. 

Fig-El shows the correlation between total p.e. and Rpsd for the neutron data. 
The nuclear recoil band just below Rpsd = 0.5 can clearly be seen, ending 
about where expected from the above calculation of maximal recoil energy. 
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At higher p.e. we also find gamma rays caused by inelastic scattering below 
Rpsd — 0.4, which is well separated from the nuclear recoils. The cluster at 
high Rpsd ~ 0.8 was caused by very sharp pulse such as Cherenkov light 
generated in the PMT window. 

4-2 Gamma-ray data 

Gamma ray data was taken with an external 137 Cs source and a lead collimator 
for the source. The trigger rate of the 137 Cs Run was 4.52 kHz (PMT1&PMT2). 
FigJT] shows the correlation between total p.e. and Rpsd for this electron 
event sample. While the average Rpsd was rather constant for nuclear recoils 
in the 252 Cf data down to low energies, it does increase for gamma-rays at 
low energies, slowly merging into the nuclear recoil region as recoil energies 
approach zero. 



5 Results 

In Fig. [6] and Fig. [7] we expect that power of Rpsd to isolate nuclear recoils 
and reject gamma-ray background is a function of the scintillation light output 
observed in our recoil event. To estimate the rejection power for gamma-ray 
background as a function of energy from our data we proceed in three steps: 
First we split our data sample into energy bands according to the total p.e. 
count observed with both, PMT1 and PMT2. We then fit Gaussian distri- 
butions to both of the Rpsd distributions, the one for nuclear recoils and 
the one for electron events, separately in each energy band. In the third step 
we estimate the contamination in the nuclear recoil sample from the tail of 
the gamma-ray Gaussian that extends beyond the mean of the nuclear recoil 
Gaussian. Using the fitted Gaussians was compared to counting events directly 
and was found to yield the same result. By integrating the tail of the Gaussian 
fitted to the electron event distribution that extends beyond the mean of the 
nuclear recoil Rpsd distribution we define the fractional electron leakage r EL 
into our nuclear recoil sample at an efficiency of 50% for nuclear recoils: 



Here G e denotes the normalized Gaussian as fitted to the Rpsd distribution for 
electron events, which has a mean \i e and a variance a e . G n , and a n stand 
for the corresponding entities derived from the fit to the Rpsd distribution of 
the nuclear recoils. 



+ inf 




(5) 
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Electron leakage was evaluated for the p.e. ranges of 100—150, 150—200, 
200-250, 250-300, and 300-400 p.e. Using a conversion factor of 20.9 p.e./keV, 
these p.e. ranges correspond to energy ranges of 4.8—7.2, 7.2—9.6, 9.6—12, 
12—14.4, and 14.4— 19.1 keV ee respectively. FigJS] shows the Rpsd distribu- 
tions for both the nuclear and the electron event samples in the energy range 
from 4.8 to 7.2keV ee . The corresponding Gaussian fits are also shown. 

FigJ9] summarizes the fractional electron leakages we measured in the energy 
ranges listed above. Horizontal bars reflect the range of electron equivalent re- 
coil energy that was used. Vertical error bars show the statistical uncertainties; 
the range between the braces has the systematic uncertainties added quadrat- 
ically. The largest contribution to the systematic uncertainty came from 1 p.e. 
determination. Open circles show the results for the study with the artificially 
reduced effective light yield of 4.6 p.e./keV ee , while solid circles show rejection 
power as a function of energy for the full measured scintillation light yield of 
20.9p.e./keV ee . 

The efficiency dependence of the nuclear recoils retention was also studied 
as shown in FigJTOl The black and red points show the rejection power in the 
energy range 4.8-7.2 keV ee and 9.6-12 keV ee for non-masked data, respectively. 
The trade off between rejection of electron events and efficiency for nuclear 
recoils retention can be evaluated from this figure. 



6 Discussion 

Our study clearly shows that PSD can be used in purely scintillation based LXe 
dark matter detectors. As expected its power depends on the effective light 
yield. Our tel measure as defined above ultimately depends on the relative 
size of two quantities: the distance of the two means fi n — fi e as compared to 
the width of the electron distribution o e . 

In FigJTTJwe compare the mean values of the Rpsd distributions for the elec- 
tron and nuclear recoil runs as a function of our energy ranges. Both the 
high yield data and the data taken with an artificially lower effective yield are 
shown, and it can be seen that the means change by less than 5% as we change 
the light yield. Yet the effect of this systematic change in \i n — /i e contributes 
only 4% to the overall efficiency loss; the larger contribution comes from the 
widening on the distribution. 

Looking at the electron distribution width a e , we tried to separate the sta- 
tistical component reflecting the p.e. statistics from an intrinsic component 
reflecting the physics of the various processes involved in generating, detect- 
ing, conditioning, and recording the signal. Fig{T2]and FigfllJ]respectivery show 
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a e and cr n , again as a function of our energy ranges. 

To estimate the statistical contribution to the width of the Rpsd Gaussians 
we used a Monte Carlo (MC) simulation. In this process we first randomly 
chose a value for p.e.tot i n the relevant range from 50 to 400 p.e.. Using the fits 
to our Rpsd means as shown in Fig{T5]we determine the proper Rpsd mean 
to use for this p.e. to t value. Using p.e.tot and its proper Rpsd, the MC obtains 
a set of p.e.prmt by sampling from the binomial distribution: 



P(p.e. prmt ) = P ,. tot C p . e . vrmt R P t6 mt (1 - RpsdY******** (6) 

From this set of p.e. prmtM c we calculate a set of R PS d,mc = P-e-prmt,Mc/P-e-tot- 
C means binominal coefficient. 

Rpsd,mc distributions are built up for both the electron and nuclear recoil 
means and fitted to Gaussians to extract fi e ,MC &e,MC, ^n,MC and a n ,MC- 
The resulting fi e ,MC and fi n ,MC are in good agreement with their measured 
counterparts \i e and fi n . 

Fig{T2] and Figff^] show the variances of a for MC and data compared for 
electron and nuclear recoil respectively. Systematically the MC data have less 
spread than the real data, but the effect is quite striking for the high light 
yield electron event data. 

In the case of electron event the waveform is dominated by recombination, 
which is slow. In case of nuclear recoil, the waveform reflects the relatively 
short lifetimes of the singlet and triplet states. Fluctuations in the energy de- 
posit along the electron track will translate into fluctuations in the recombina- 
tion time scale and might be one mechanism by which additional fluctuations 
can propagate into our Rpsd parameter. In an effort to quantify the extra 
contribution Oi n t r insic to the width of our distribution that is evident in our 
data we calculated: 



'intrinsic — y^e ® e,MC {J J 

The results are summarized in Fig. [15j The observed consistency in all over 
energy range gives support to the hypothesis that we are seeing the effects of 
processes at the light emission stage, as that part was not changed; our light 
yield was only effectively reduced by a mask. Intrinsically the light yield in 
the low and high yield data runs was still the same. 
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7 Conclusion 



PSD in liquid xenon was studied at energies relevant to dark matter searches. 
A significant difference in pulse shape between nuclear recoil and electron 
events was exploited in a high light yield setup in which 5keV ee energy deposit 
produce 100 p. e. in PMT signal output. 

At high light yield (20.9 p.e./keV) our Rpsd parameter allows a rejection of 
electron events to a level of 7.7±l.l(stat)±g g(sys) x 10~ 2 in the energy range 
between 4.8-7.2 keV ee with a 50% efficiency to retain nuclear recoil events. In 
the energy range 14.4-19.1 keV ee electron events were reduced by more than 3 
orders of magnitude with the same efficiency for nuclear recoils. 

The dependence of this rejection power on photon statistics was also stud- 
ied. At low effective light yield (4.6p.e./keV), a rejection of electron events 
to 2.4±0.2(stat)±°J( S ys)xlO- 1 was demonstrated for the energy range 4.8- 
7.2 keV ee with a 50% efficiency for nuclear recoil events. 

In our MC replication of the experimental results the width of the Rpsd 
distribution for nuclear recoils is almost exhausted by the expected purely 
statistical contribution to that width, while that is clearly not the case for 
electron events. This excess width was tentatively interpreted as stemming 
from fluctuations inherent to the energy deposit of electron events. 
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Fig. 1. A schematic view of the detector used in the PSD measurement. 
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Fig. 2. A schematic diagram of the data acquisition system. 
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Fig. 4. A typical p.e. distribution of Co data for PMT2. The light yield was reduced 
to 4.6p.e./keV ee by a mask. 
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Fig. 5. TOF distribution of 252 Cf Run. 
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Fig. 6. The scatter plot between total p.e. and Rpsd for TOF from 15 to 30nsec 
in the 252 Cf data. 
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Fig. 7. The 137 Cs scatter plot between total p.e. and Rpsd- 
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Fig. 8. The Rpsd distribution in the energy range between 4.8-7.2 keV ee . The red 
and black lines show the 137 Cs data and the 252 Cf data, respectively. 
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Fig. 9. The rejection power distribution. The open and filled circle points show the 
masked(4.6 p.e./keV) and non-masked (20.9 p.e./keV) data. Error bars shows the 
statistical error. Braces show the error including systematic uncertainty added in 
quadrature. 
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Fig. 10. The efficiency dependence of nuclear recoils retention for rejection power. 
The black and red points show the rejection power in the energy range 4.8-7.2 keV ee 
and 9.6-12 keV ee for non- masked data, respectively. 
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Fig. 11. Mean value of the Rpsd as a function of observed energy in unit of keV ee 
for 137 Cs runs and 252 Cf runs. The filled and open triangle points show the values 
of 137 Cs runs for the non-masked PMT1+PMT2 (20.9 p.e./keV) and for the masked 
PMT2 (4.6 p.e./keV), respectively. The filled and open circle points show the values 
of 252 Cf data for the non-masked PMT1+PMT2 (20.9 p.e./keV) and for the masked 
PMT2 (4.6p.e./keV), respectively. 
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Fig. 12. Sigma of the Rpsd distribution as a function of observed energy in unit 
of keV ee for 137 Cs data. The filled circle and filled triangle points show the values 
for the non-masked PMT1+PMT2 data (20.9 p.e./keV) and for the masked PMT2 
data (4.6 p.e./keV), respectively. The open circle and open triangle points show the 
MC simulation of the non-masked PMT1+PMT2 (20.9 p.e./keV) and the masked 
PMT2 (light yield was reduced to 4.6 p.e./keV), respectively. 
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Fig. 13. Sigma of the Rpsd distribution as a function of observed energy in unit 
of keV ee for 252 Cf data. The filled circle and filled triangle points show the val- 
ues for the non-masked PMT1+PMT2 (20.9 p.e./keV) and for the masked PMT2 
(4.6 p.e./keV), respectively. The open circle and open triangle points shows the MC 
simulation of the non-masked PMT1+PMT2 (20.9 p.e./keV) and the masked PMT2 
(light yield was 4.6 p.e./keV). 
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Fig. 14. The relation between Rpsd and total p.e. The black and red lines show the 
fit function of 252 Cf and 137 Cs. The error bar of Cs represent the 1 sigma from 
the mean value. 
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Fig. 15. The (^intrinsic distribution of electron events. The open and filled circle 
points show the values of the masked and the non-masked data, respectively. 
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